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We investigate the stochastic dynamics of a particle in the presence of a modulated sinusoidal 
potential. Using the time derivative of the winding number, we quantify the particle's motion 
according to its running time, the time it runs monotonically to the left or right. For a range of 
model parameters, we show that, in the overdamped regime, the particle's motion in this modulated 
washboard potential exhibits stochastic resonance. We briefly suggest possible applications of our 
results, including the amplification of signals for measurement devices and in stimulated tunneling 
of Bose- Einstein condensates. 

PACS numbers: 



INTRODUCTION 

Stochastic resonance (SR), originally proposed in ref- 
erence to two equally-favored climate states, has only 
more recently been discussed as a more general concept 
independent of an explicitly bistable potential function 
(see (H for a review). It has now been identified in 
a variety of different contexts, from lasers and optical 
traps to neurological systems and biochemical reactions. 
Quantum versions of stochastic resonance have also re- 
ceived attention, especially in superconducting quantum 
interference devices (SQUIDs) and nanomechanical res- 
onators due to their location near the boundary of clas- 
sical and quantum dynamics. It has been shown that SR 
can be an efficient amplification mechanism, and similar 
ideas are of interest to controlled diffusion, for instance 
by applying various noise models in symmetric or asym- 
metric periodic potentials to obtain directed transport 
(such as described by "Brownian motors" 0). 

The interest in nonlinear resonance phenomena in pe- 
riodic potentials is not due solely to models of diffusion 
and transport, as such potentials underlie other inter- 
esting physical systems, for example, sphalerons in high 
energy physics, which have been invoked in the context of 
the electroweak phase transition Q . Sticking to the con- 
crete language of diffusion, many authors have discussed 
SR in the basic sinusoidal potential, starting with the 
simplest (static, symmetric) case 0, IH , implementing a 
static bias [l, 0, B, finally tackling actively modu- 
lated systems [ol, iiol . [HI . Ilij . Borromeo and Marchesoni 
recently generalized stochastic resonance to dynamically- 
defined states in a tilted washboard potential . 

One case not yet treated in the literature is a rocking 
washboard, in which a periodic potential is subjected to 
a global modulation, producing alternately positive and 
negative biases. This system offers a simple model for dif- 
fusion in the presence of agitation, such as experienced 
by adatoms or particulates on a vibrated surface, or even 
electrons on the surface of a dielectric and subject to 
a radiation field. Periodic potentials are also known in 
particle physics, so the possibility of an oscillatory bias, 
e.g., through a background scalar field Q, further rec- 



ommends the relevance of this model. In the next section 
we introduce the model and explore some of its proper- 
ties. In section 3 we discuss our results, while in section 
4 we present a summary. 



MODEL 

The equation of motion is the usual one associated with 
studies of Brownian motion and stochastic resonance. 



x = (,{t)--fx-V'{x,t) 



(1) 



in rescaled units. The first two terms on the right 
represent the thermal bath at temperature T via the 
fluctuation-dissipation theorem: Gaussian white noise 
of zero mean (^(t)) = and 5 correlation function 
{£,{t)^{t')) = 2jkTS{t — t') combined with viscous damp- 
ing 7. k is Boltzmann's constant. The last term is the 
potential, which has both a static and a periodic driving 
component 

V{x,t) = -(A/27r)cos(27rx) -^xFcosnt. (2) 

Previous work [13, llH has maintained a constant 
bias throughout any time-dependent modulations of bar- 
rier heights, so that the driving component appeared as 
F{t) = Fo + FicosQt for Fq > 0. Our study departs 
by setting Fq = 0, thereby tilting the whole washboard 
rather than merely modulating its corrugation. 

We consider a low frequency input and the overdamped 
limit, in which the equation of motion ([T]) can be adia- 
batically reduced to 



-A sin(27rx) + F cos nt + ^{t) 



(3) 



with t normalized to absorb 7. The potential was defined 
so that integer values of the position correspond to the 
winding number co. 

In the limit of slow modulation, 1] <C 1, transitions 
between wells are thermally induced for any value of the 
damping. At various points in its evolution, the potential 
has a nearly constant bias in the positive or negative 
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direction, and the difference induced in the right and left 
barrier heights is given by 

V{l/2,t) - y(-l/2,t) = Fcosnt. (4) 

For phases near Q = 7r/2,37r/2, regular diffusion across 
an unbiased periodic potential applies, while around 
1] = 0, TT, the dynamics are effectively those of activated 
diffusion. In the latter case, the extensive characteri- 
zation of the jump length and time statistics given by 
[l, [lH for the tilted washboard potential may be straight- 
forwardly applied. Jump lengths and times, respectively 
defined as the number of wells and the total period of 
monodirectional motion, are controlled by the ratio F/7, 
as discussed in [l3| . 

The activated behavior near maximum bias will only 
dominate the overall statistics of the system if the particle 
is effectively locked into a specific well for the remainder 
of the phase. Letting E~ = A/tt — F/2, the reduced 
barrier height in the direction of the bias, left at 1] = 
and right at 1] = tt, we may derive an intuitive con- 
straint in the spatially-diffusion-controlled regime (large 
damping in our normalization) via Kramers rates [Iq . 
On the one hand, the particle should experience the 
tilted washboard potential for longer than the proba- 
ble escape period, Tmin < Tq/A with Tmin = l/vk = 
{27r J /ujquJ})) ex.p{PE~) the inverse of the Kramers rate 
near maximum bias and the quarter period being a rea- 
sonable choice, (c^o and 00^ correspond to the oscillating 
frequencies at the metastable well and top of the bar- 
rier, respectively, and (3 = 1/kT.) On the other hand, 
the particle should also have a suppressed probability 
of escape near zero bias, and hence T^ax > Tq/A for 
Tmax = {27rj/uJoUJ})) ex.-p{(3A/7r). Rearranging and com- 
bining the inequalities, we find 

< {kT/A) ln(4ao) + I/tt < F/2A (5) 

where ao = CIj/coq is a dimensionless time parameter in 
which we used the fact that = uoq for the sinusoidal 
potential. Because this inequality locates the region of 
parameter space in which we expect the activated behav- 
ior to dominate the overall dynamics, we clearly desire 
F/A as large as possible. 

A very similar constraint, more accurate in the case of 
low damping, can be derived from the Kramers rate in the 
energy-diffusion-limited regime, ^I[E) <C kT (in which 
I{E) is the energy-dependent fiux in action space). An 
extra step is required in combining the inequalities due 
to the different values I{E~) and /(A/tt), 

HAM^e^^^- < 1, (6a) 

I^E-)^e^^/- > e^^l\ (6b) 

We then use the monotonicity of I{E\ from which 
1{E-) > /(A/tt), to write 

< (kT/A) \n{4au) + I/tt < F/2A. (7) 



We retain the larger fiux I{E~) in order to strengthen 
former inequality (|6a|). since it expresses the suppression 
of the fiux during the unbiased parts of the period. For 
the effect to become pronounced, the inequality should 
clearly be as strong as possible. However, in the context 
of stochastic resonance, small input biases are of greater 
interest, and in the following F < A. 

In a biased periodic potential, overdamped systems 
represented by (|3]) can exhibit a consistently running 
state only for F — 2^4, as the barriers are completely 
suppressed in this case. However, many results for diffu- 
sion in a tilted periodic potential are quite general [13], 
and directed diffusion still arises with an average velocity 
given by 

{kTh){l-eM-F/kT) 
/ Si dxdyeMiy{x) - V{x - y) - yF)/kT] ' 

Considering shorter time dynamics, the modulation of 
the barrier heights introduces a phase bias for diffusive 
transitions. For small F, a simple correlation should be 
apparent in the residence time distribution, as jumps be- 
tween wells will tend to occur in phase with the modula- 
tion. 

We therefore expect a strong correlation between the 
direction of transport and the phase of the forcing. Given 
that transitions are suppressed when the biasing force 
F{t) = F cosQt is near zero, transitions against the prof- 
fered bias are even more strongly suppressed (though ve- 
locity inversions occur less often than initially expected 
in a statically biased potential according to [8|). 

We observed a dynamic resonance between the time 
derivative of the winding number w and the modula- 
tion F{t), now understood as an input signal. Defining 
"left-running" and "right-running" states corresponding 
to It; < and w > 0^ respectively, the formalism and lan- 
guage of stochastic resonance can be applied. All of the 
generally accepted ingredients are present: (i) an activa- 
tion barrier, E~ = 2 A — F/2; (ii) a (weak) input signal, 
the global modulation of the potential, F{t) in and 
(iii) a source of noise, namely the stochastic term in ([1]). 

RESULTS 

The most direct connection to previous observations 
of SR lies in residence times, defined as the time spent 
in a given well: we found that the distribution of res- 
idence times for the left- and right-running states — or, 
more appropriately to our situation, the running times — 
are peaked around the half integer multiples of the forcing 
period, (n + ^)Tq. Running times for various values of 
the temperature are plotted in figure (TJ and the qualita- 
tive similarity to plots for resident times such as found 
in [1] is easily verified. 

We note that the usual resonance in the residence time 
distribution (time spent localized in a given well) can still 
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FIG. 2: The Fourier transform of position x in an over- 
damped system, driven at Q — 10~^ . The primary peak at 
the driving frequency is clear, accompanied by an extra peak 
at the first harmonic, SQ, exactly as predicted by Pa- 
rameters: A = 1, F = .3, 27vkT = .325 



FIG. 1: The residence (top) and running (bottom) time dis- 
tributions in the overdamped model for three temperatures 
27TkT/A = 0.25, 0.35, 0.45. The vertical axis denotes the num- 
ber of instances of a measured period, normalized on the hori- 
zontal axis by the forcing period, Tq. The residence time also 
displays a (barely-discernable) resonance, showing a prefer- 
ence for transitions on both integer and half-integer multiples 
of the period at sufficiently low temperature. The distribution 
of the jump length Aw in the inset shows the expected expo- 
nential die off [5]. Parameters: A = lO'^, F = 2000, Q = 10 



be identified at low temperatures, although not as clearly 
due to the scale of the figure. Below, we will present 
results in the underdamped regime, where the resident 
time resonance is more apparent (see Fig. 4). 

With 7 normalized to 1, F/7 is always less than 
(the threshold for the locked regime so jump lengths 
display the expected exponential distribution, as can be 
seen in the inset of figure [T] [5i] . 

Note that 'running' refers only to the monotonicity of 
the winding, as the particle almost always resides a time 
comparable to the relaxation time cJq^ even when only 
transiting a well. For this reason, in order to obtain the 
clean results displayed in figure [U we considered a change 
in winding number to occur only after the particle begins 
its descent into the adjacent well, not merely when it has 
surmounted the barrier. This definition is the cause of 
the low frequency cutoff apparent upon close inspection. 

Further support for characterizing the running time 
as a resonance comes from the Fourier transform of x. 
The signal-to-noise ratio should show a maximum in the 
temperature range suggested by ([5j) . Figure [2] displays 
the Fourier transform of a system in the suggested range 
and furthermore exhibits the first harmonic at as pre- 
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FIG. 3: The normalized diffusion coefficient D/Dq plotted 
against temperature, obtained from the regression of (x^) 
from 1000 particles. Parameters: A = 1,F = .3,Tq = 1000 

dieted by 0. 

In [l2[, it is noted that diffusion is somewhat sup- 
pressed in a tilted periodic potential (as compared to 
Einstein diffusion I^o = kT/j)^ and the global modula- 
tion may be expected to suppress diffusion even further 
by introducing an extra coherence in the transport of 
particles in the potential. As can be seen in figure [3] we 
have investigated this fact, finding a very significant sup- 
pression, which decreases with increasing temperature. 

In underdamped systems 7 <C cjq, the left- and right- 
running states have an obvious counterpart in the jump 
times and lengths discussed in However, the corre- 
lation of w to the phase of the drive cannot be properly 
called a resonance phenomenon, as inertial effects rein- 
force rather than suppress the peaked distribution, as 
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FIG. 4: Residence (bottom) and running (top) time distri- 
butions for an underdamped system at temperatures 27vkT = 
0.2, 0.25, 0.3. The response of w does not show a resonance- 
Uke behavior, as the distribution remains sharply peaked at 
the half integer multiples even while temperature increases. 
The residence time is also dramatically peaked near zero, 
though the higher harmonics here also gain power as temper- 
ature increases. The steadiness of the distribution, however, 
highlights that it is entirely an inertial effect. Parameters: 
A = l,F = 0.3,Tn = 100 

demonstrated in figure [H Physically, if 7 is of the same 
order as 2^7, sufficient kinetic energy may be conserved 
through the half period of opposing bias such that no ve- 
locity inversion occurs (8| . Evidence for this explanation 
is derived from the Fourier transform of x: the power 
in peaks at higher harmonics should generally increase 
with temperature, as lower frequency velocity inversions 
become more common. Simulations support this, to a 
point, as we find that the peaks actually become satu- 
rated beyond a certain temperature, as shown in figure [3 
Although the saturation indicates a decrease in signal- 
to-noise ratio, it does not produce a profile sufficiently 
peaked as to justify its characterization as a resonance. 



SUMMARY 

We described a dynamic resonance exhibited by an 
overdamped particle in a rocking washboard potential. 
Defining "running" states according to the sign of 
whereby a particle runs monotonically to the left or right 
before reversing its direction, we were able to characterize 
the phenomenon as a new kind of stochastic resonance. 
Numerical results are consistent with the range of pa- 
rameters derived from naive expectations. We have also 
shown that no such resonances are possible in the under- 
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FIG. 5: The power in the delta peaks of the first through 
fourth harmonics, (2n + 1)^ for n = 1..4. The power in 
each peak increases nonlinear ly with temperature, indicating 
an increase in signal to noise ratio — behavior contrary to the 
resonant character of SR. Each harmonic appears to saturate 
around 2iTkT = 0.4. The inset displays the Fourier transform 
of an experimental run at 27ikT = 0.35 in which the peaks at 
the higher harmonics can be clearly identified. Parameters: 
A = 1, F = 0.3, 7 = 0.015, Tn = 100 

damped regime. 

Looking ahead, our results are of interest to measure- 
ment considerations in the context of bifurcation ampli- 
fication [1^ , whereby a small input signal is enhanced so 
as to effect the distinction between the (bistable) dynam- 
ical states arising just beyond the bifurcation threshold 
of a nonlinear system. This immediately suggests a pos- 
sible direction in which these ideas may be taken: instead 
of F{t) being an unvarying input mode, the particle may 
be coupled to a second system whose oscillations provide 
the input signal. From this perspective, ([T]) represents 
the high mass-ratio limit of this second system, so that 
the observed variable (modeled here) effects no significant 
back-reaction on the input signal. For example, starting 
from the equations of motion for a bilinearly coupled pair 

mixl = -Vl{xi) - 71X1 + Xx2 
m2xl = -V2{x2) - 72^2 + Axi 

we may consider the perturbative case in which the so- 
lution X2{t) is harmonic, being limited to a region near 
one of its potential minima. Secular perturbation the- 
ory on the full system produces the expected condition 
for the validity of this Ansatz, namely A/m2 <C 1. Then 
with ^ a model of such a bipartite system, the dynami- 
cal resonance described above becomes a large amplitude 
signal for the oscillation of the underlying system {X2) be- 
tween two of its states. These states may be dynamical 
or static, and the coupling engineered accordingly. 

A second connection can be made to recent advances in 
quantum tunneling dynamics. Experiments were recently 
undertaken in which resonantly enhanced tunneling was 
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clearly observed in a Bose-Einstein condensate in a tilted 
washboard potential |20|] . Although the proper quantum 
mechanical calculations should be carried out in detail, 
the classical results here suggest that by tuning the oscil- 
lation of the washboard to the tunneling rate, a second 
order resonance should arise. The loss of particles to the 
continuum following the first tunneling event could be 
drastically reduced by reversing the tilt on the potential. 
The discretization of the energy levels provides an inher- 
ent suppression of tunneling events during the inversion 
of the bias. Testing this hypothesis might more easily 
be achieved in a circularly symmetric optical potential 
such as was used in [21], so that the modulated bias is 
alternately positively and negatively radial. 

This paper demonstrates the utility of stochastic reso- 
nance as a general model of phenomena for which obvious 
localized states are not necessarily presented by the form 
of the potential. Clearly, with ideas progressing in such 
directions as bifurcation amplification [19] and the abil- 
ity to easily construct and manipulate periodic potentials 
thanks to laser trapping, the language and the theory to 
describe the possible dynamics should develop as well. 
Stochastic resonance in the context of more generalized 
states will almost certainly appear and be fruitfully dis- 
cussed in ways similar to those presented above. 
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